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We perform Monte Carlo simulations of the Ising spin glass at low temperature in three dimensions 
with a ± J distribution of couplings. Our results display crossover scaling between T = behavior, 
where the order parameter distribution P(q) becomes trivial for L — > 00, and finite-T behavior, where 
the non-trivial part of P(q) has a much weaker dependence on L, and is possibly size independent. 
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I. INTRODUCTION 



Several papersoQ have recently studied the Ising spin 
glass in three dimensions with a Gaussian distribution 
of bonds at low and zero temperature. From data ob- 
tained on small sizes these papers deduce that the order 
parameter distribution function, P(q), is non-trivial at 
finite-T, i.e. in addition to two peaks, symmetric about 
q = 0, there is also a continuous part betwjeen the peaks 
whose weight does not decrease with sizecl. This indi- 
cates the existence of a nontrivial energy landscape, i.e. 
of macroscopic excitations, involving a finite fraction of 
the system, that cost a finite energy in the thermody- 
namic limit. This aspect of the results is consistent with 
the replica symmetry breaking picture of Parisij. By 
contrast, the droplet theory&El predicts that the weight 
in the continuous part of the distribution should vanish 
like L as the (linear) size of the system L increases, 
where 9 is a positive exponent. In both theories, because 
the ground state is unique (apart from inverting all the 
spins), it follows that the weight in the "tail" of the dis- 
tribution tends to zero (proportional to T) as T — > and 
the positions of the peaks tend to ±1. The purpose of 
this paper is to see how these results are modified for 
a spin glass with a bimodal distribution (also called the 
±J distribution), where the interactions have values ±1, 
where there is a large ground state degeneracy and a fi- 
nite ground state entropy per spin. 

One might possibly imagine that, since the system with 
the ± J distribution has a finite ground state entropy, its 
behavior at zero temperature would be similar to that of 
a model with continuous distribution at finite-T. If 
were true then, according to the numerical results,! 
P{q) would be non-trivial at T = whereas according to 
the droplet theory P(q) would be trivial. 

However, -this notion has been contested by Krzakala 
and MartinO (referred to henceforth as KM) who argue 
that entropy effects cause one "valley" in the T = en- 
ergy landscape of the ± J model to dominate and conse- 
quently the weight in the tail vanishes like T~ A , where A 
is a positive exponent (discussed below), even if the en- 
ergy landscape is non-trivial. At finite-T, KM argue that 
the weight is finite for large L, so, by implication, there 



must be a crossover at some scale L C (T) from the L~ x 
behavior for L < L C (T) to a value independent of L at 
larger sizes. One can also generalize the KM argument to 
the droplet model, in which case there is still a crossover, 
between L~^ XJr9 ^ behavior at smaller L and L~ e behavior 
at larger L. Overall, in the KM scenario, the only dif- 
ference between the continuous and the ± J distributions 
for L 3> L C {T) is that the position of the peaks in P{q) 
are different for T — > 0. Denoting the peak positions by 
±go, then one has go < 1 for the ± J distribution whereas 
qo = 1 for a continuous distribution. 

Here we display, we believe for the first time, the 
crossover between T = and finite-T behaviors. Further 
motivation for our work is to clarify conflicting results 
for ground state properties. Berg et al£j used a mul- 
ticanonical Monte Carlo technique to determine P(q) at 
T = finding results consistent with trivial behavior with 
A = 0.72 ± 0.12 (but also not ruling out the possibility 
of nontrivial behavior). HartmannEJ used a genetic opti- 
mization algorithm finding initially a nontrivial P(q), but 
the results were biasecfij because the degenerate ground 
states were not sampled with equal probability. Subse- 
quently Hartmannllia developed an improved method and 
found a trivial P(q) with A = 1.25 ± 0.05, and suggested 
that this supports the ^droplet picture. Very recently 
Hatano and GubernatisE-3 (referred to as HG) have per- 
formed a "bi-variate multi-canonical" Monte Carlo study, 
finding that T(0) drops dramatically at low-T as L in- 
creases. Though they do not extract the exponent A, 
from the figures in their paper, it appears that A is sig- 
nificantly larger than Hartmann's value. They too argue 
that their results provide evidence for the droplet pic- 
ture. However, Marinari et al.E_3 have recently claimed, 
on the basis of their own simulations, that the results of 
HG are not equilibrated and thei^conclusions are there- 
fore invalid. Finally, recent workt^l finds a nontrivial en- 
ergy landscape and also, apparently, a nontrivial P(q) at 
T = 0. It therefore seems useful to try to decide between 
these different results. Our data at the lowest tempera- 
tures imply a trivial P(q) at T = and .our estimate for 
A is consistent with that of Berg et al.EJ but not with 
that of HartmannEj or HG. 

The Hamiltonian is given by 
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where the sites i lie on a simple cubic lattice in dimen- 
sion d = 3 with N = L 3 sites (L < 10), Si = ±1, and 
the Jjj are nearest- neighbor interactions taking values ±1 
with equal probability. We do not apply the constraint 
Y^Uj) Jij = 0, which is imposed in some related work. 
However, we expect that the crossover from T = to 
finite-T behavior will be similar in the two models. Pe- 
riodic boundary conditions are applied. We focus on the 
distribution of the spin overlap, q, where 
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in which "(1)" and "(2)" refer to two independent copies 
(replicas) of the system with identical bonds. 
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FIG. 1. An equilibration plot for L = 8,T = 0.20, for 
the second and fourth moment of P(q), and for x(l/2), the 
average of P(q) over the interval \q\ < 1/2. For better viewing, 
the data for (q 4 ) and x(l/2) have been shifted upwards by 0.14 
and 0.67, respectively. For each value of N aweep , the averages 
were measured /3 MC sweeps. 

Simulations of spin glasses at very low temperatures 
are now possible, at least for modest .sizes, using the 
parallel tempering Monte Carlo methodlJO, where one 
simulates replicas of the system at Nt different temper- 
atures. Here, we need two copies of the system at each 
temperature to calculate q, so we actually run 2 sets of 
Nt replicas. Mle also gain a large speed-up by using mul- 
tispin codingij to store each spin or bond as a single bit 
rather than a whole word. 

In earlier workEl for the Gaussian distribution we were 
able to use a special relationship between certain vari- 
ables to check for equilibration, but this is not applicable 
here. We therefore investigate whether various quantities 
have become independent of simulation time when plot- 
ted on a logarithmic scale. Fig. [j] shows an example for 
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TABLE I. Parameters of the simulations, 
number of samples (i.e. sets of bonds), N svacp is the total 
number of sweeps simulated for each of the 2Nt replicas for a 
single sample, Nt is the number of temperatures used in the 
parallel tempering method, and T m i„ is the lowest tempera- 



ture simulated. N si 



=6336 for L = 8 and T > 0.35. 



L = 8,T = 0.20 indicating that the data seems to have 
saturated. 

In Table [|, we show the simulation parameters. The 
lowest, temperature simulated, T m i n , has to be compared 
witbEJ T c w 1.15. For each size the largest temperature is 
2.0. The set of temperatures is determined by requiring 
that the acceptance ratio for global moves is 0.3 or larger. 
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FIG. 2. Data for the overlap distribution P(q) at T = 0.20. 
The vertical scale is logarithmic to better make visible the 
peak at large q and the tail down to q — 0. We only display 
some of the data points as symbols, for clarity, but the lines 
connect all the data points. This accounts for the curvature 
between neighboring symbols. 

Figs. [2] and || show data for P(q) for different sizes at 
T = 0.2 and T = 0.35. One can see that the weight in 
the tail tends to decrease initially with increasing L, espe- 
cially at lower T, but for T = 0.35 the data seems to sat- 
urate at larger L. For T = 0.5 (not displayed) the weight 
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in the tail saturates already at L — 4. This can be seen 
more clearly in Fig. [|, which shows x{\/2) as a function of 
L for different temperatures, where x(q) = /' P(q')dq' 
so x(l/2) is the average of P(q) from —1/2 to 1/2. We 
give data for x(l/2) rather than P(0) because the statis- 
tics are better and also so we can compare directly with 
other work. For L = 4 and L = 6, the data at T = 0.05, 
not showed in Fig. [|, are superimposed to the data at 
T = 0.2, indicating that we have reached the true T = 
behavior. For L = 8, the data at T = 0.2 may be one 
or two standard deviations larger than the T = value. 
Furthermore, the average energgi-at T = 0.2 agrees with 
the ground state results by PalEiL From a power law fit 
of the data in Fig. || at T = 0.2 we estimate 

A = 0.9 ±0.1. (3) 




FIG. 3. Same as for Fig. | but at T = 0.35. 

Generalizing the KM argument to a scenario described 
by an exponent 6, we expect that at finite-T there will be 
a crossover between the L~( A+e ) behavior for L smaller 
than some length scale L C (T), and the L~ e behavior (or, 
for 8 = 0, an L-independent value proportional to T), 
at scales larger than L C (T). In the more general case, 
assuming scaling one has L C (T) ~ T _1 / A and 

a;(l/2) =TL- 6 f{LT l ' x ), (4) 

where / is a scaling function. 

A scaling plot appropriate to this behavior, for 9 = 
and A = 0.9, is shown in Fig. |[ where one can see that 
the data collapse fairly well. The data in Fig. || increase 
with increasing L for T > 0.8, due to the vicinity of 
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FIG. 4. Log-log plot of :r(l/2), the average of P(q) over 
\q\ < 1/2, against L. 

T c , where x(l/2\ ~ L p / U andEl (3/v ~ 0.3. One may 
therefore argueE3 that the observed saturation between 
T = 0.35 and T = 0.65 is a finite size effect and that at 
larger sizes there will be a second crossover to the L~ e 
behavior. We cannot exclude this possibility, though we 
note that T — 0.35 is quite far from T c and that a scaling 
plot as in Fifr,|but with 9 = 0.2 is significantly worse. 

Hartmannta computed ar(l/2) as a function of L at 
zero temperature and found that a power law fits well 
the data with an exponent A = 1.25 ± 0.05, which dis- 
agrees with our estimate. Ou*-,value for A does, however, 
agree with that of Berg et alJlS who find A = 0.72 ± 0.12. 
In addition, our raw data for x(l/2) is consistent, with 
(though more accurate than) that of [Berg et al.E2l, but 
is inconsistent with that of Hartmannllj for L > 4. For 
example, for L = 6 we find at(l/2) = 0.095 ± 0.002, while 
Hartmann finds x(l/2) — 0.083±0.005. We note however 
that Hartmann's method, unlike (properly equilibrated) 
Monte Carlo simulations, is not guaranteed to sample all 
the ground states with equal probability. 

Our results for P(q) at low-T are also in marked dis- 
agreement with HG. For example, HG report a P(q) 
which is lower than 0.03 in the interval \q\ < 0.1 for L = 8 
and T = 0.3, while our average of P(q) over this interval 
is between 0.066 ± 0.004 (our value at T = 0.275) and 
0.081 ± 0.004 (our value T = 0.35). HG observe a pro- 
nounced decrease of P(q) with L even at T = 0.5, where 
our data clearly saturate. We also computed the Binder 
cumulant, which agrees with Ref. |l5] but disagrees with 
HG. This suggests that the simulations of HG are not 
correctly equilibrated, as discussed in detail in Ref. IHjL 
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FIG. 5. The scaling behavior of x(l/2) expected from 
Eq. (|) with = 0. For L » L C (T) ~ T~ 1/A , x(l/2) is 
independent of size, while for L <C L C (T), x(l/2) varies as 
The dashed line has a slope of —0.9. 



KM give arguments that A should equal d s /2 where d s 
is the fractal dimension of the surface of the large-scale 
low-energy excitations which give rise to a nontrivial en- 
ergy landscape. However, one expects that d s > d — 1 
which is barely satisfied by the estimate in Eq. (^) which 
corresponds to d s — 1.8±0.2. Furthermore, for the Gaus- 
sian distribution, d s is significantly larger than this value. 
For example, Ref. | finds d s = 2.58 ± 0.02. While it is 
possible that d s could be different for the Gaussian and 
±J models, our results suggest that A ^ d s /2, and that 
there may be corrections to the argument of KM. 

To conclude, results from simulations on small sizes 
indicate that the order parameter distribution of the ± J 
Ising spin glass is trivial at T — but, at least for quite 
small sizes, is nontrivial at finite-T in agreement with the 
conclusions of KM. We have also demonstrated crossover 
scaling between the zero-T and finite-T behaviors. We 
expect similar results in other models with a discrete dis- 
order distribution, and indeed this is what we find in pre- 
liminary unpublished data for the ± J Ising spin glass in 
d = 4. Whether these conclusions are still valid in the 
thermodynamic limit remains an open question. How- 
ever, we emphasize, quite generally, that a trivial P(q) 
at T — does not, in itself, imply evidence for the droplet 
model since this is also expected if P(q) is nontrivial at 
finite-T, as pointed out by KM. 

After this work was submitted we received a paper by 
Hed et alB3, in which, based on a different analysis from 
ours, they claim that P(q) is non-trivial at T = 0. 
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